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Abstract 

The  Born-von  Karman  model  for  a  chain  polymer  crystal 
previously  studied  by  Stockmayer  and  Hecht  is  re-examined 
using  only  analytic  approximations  rather  than  numerical  methods. 
The  analytic  approach  brings  out  many  peculiar  properties  of 
the  chain  polymer  model  that  did  not  appear  in  the  numerical 
treatment. 

One  branch  of  the  frequency  distribution  g(v)  is  shown 
to  be  proportional  to  v£  for  the  smallest  values  of  v,  approxi¬ 
mately  proportional  to  v^2  for  slightly  larger  v,  to  for 

still  larger  v  and  to  v~1//2  in  still  a  fourth  range  of  small 
values  of  v.  An  accurate  graph  of  g(v)  is  constructed  for  the 
entire  range  of  frequencies  using  values  of  the  force  constants 
suggested  by  Stockmayer  and  Hecht,  g(v)  is  shown  to  have 
approximate  singularities  of  a  type  not  anticipated  by  van  Hove 
in  his  broad  treatment  of  singularities  for  general  systems. 

This  anomalous  behavior  results  from  having  strong  valence 
forces  resisting  the  bending  of  bond  angles,  A  classification 
and  description  of  various  kinds  of  singularities  that  may 
arise  for  systems  with  strong  valence  forces  are  given  in  an 


appendix. 
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1*  lQtjrpdU9t;pfl 

The  specific  heats  of  lamellar  crystals  have  been  the 

subject  of  considerable  study,  both  experimental1*  and  theoret- 
2 

leal  *,  in  the  last  few  years.  Chain  polymer  or  fibrous 
crystals,  on  the  other  hand,  have  received  much  less  attention 
despite  the  fact  that  they  should  show  anomalies  at  least  as 
interesting  as  those  of  the  lamellar  crystals. 

Stockmayer  and  Hecht^*  considered  a  very  simple  rectangu¬ 
lar  lattice  as  a  model  for  a  hypothetical  chain  polymer  system 
and  analyzed  the  spectrum  using  the  numerical  scheme  of 
Blackman.  The  model  which  they  proposed  is  sufficiently  simple, 
however,  that  much  of  the  analysis  can  be  done  analytically 
instead  of  numerically  with  the  result  that  certain  features 
of  the  spectrum  and  of  the  specific  heat  that  were  hidden  in 
the  numerical  treatment  now  become  more  apparent.  It  is  the 
purpose  of  the  following  to  extend  the  analysis  of  this  model 
and  to  point  out  which  of  these  features  are  representative, 
at  least  qualitatively,  of  a  whole  class  of  chain  polymer 
orystals  and  which  features  are  peculiar  to  this  particular 
model* 

We  shall  use  the  same  basic  notation  here  as  in  S-H  with 
the  following  exception;  since  we  shall  be  using  at  all  times 
the  reduced  force  constants  and  frequencies  denoted  in  S-H  by 
subscripts  r,  it  is  convenient  to  suppress  this  index.  The 
a,  0,  x,  v  etc.  used  here  shall  be  interpreted  as  the  ar,  Pr, 
xr,  vr  etc.  of  S-H.  The  reader  is  referred  to  this  paper  not 


Nonr  562-08 


3 


only  for  the  notation  hut  for  the  discussion  of  the  model,  and 
the  general  formulation  of  the  problem.  We  would  like  to  point 
out,  however,  certain  significant  features  in  this  paper  that 
might  otherwise  pass  unnoticed, 

S-H  apparently  were  the  first  to  incorporate  into  the 
Born-von  Karman  model  a  force  constant  for  bending  of  valence 
angles  even  though  it  has  long  been  common  practice  to  include 
such  forces  in  the  analysis  of  the  vibrations  of  complex  mole¬ 
cules.  The  existence  of  the  force  constant  x  particularly  with 
x  »  a,  y  Is  considerable  consequence  as  will  be  shown  below. 
The  combination  of  rectangular  geometry  and  central 
forces  only  between  chains  leads  to  the  result  that  in  Eq,  (2,8) 

of  S-H.  v2  (v2)  depends  only  very  weakly  upon  cp  for  y  «  1* 
l  3r  2 

This  situation  also  has  a  pronounced  effect  upon  the  spectrum 
and  is  in  some  respects  similar  to  the  effect  seen  in  the 

L 

analysis  of  the  simple  cubic  crystals  !  in  which  second  nearest 

neighbor  forces  are  necessary  to  prevent  the  spectrum  from 

becoming  essentially  one  dimensional  in  form.  Any  consequences 

resulting  from  y  «  a  are  closely  tied  to  the  rectangular 

geometry  of  the  crystal  and  are  not  to  be  considered  as  typical 

of  all  chain  polymer  crystals  whereas  the  consequences  of 

x  »a,y  probably  are  typical  as  are  also  those  of  »  y  , 

The  peculiar  effects  of  the  rectangular  geometry  arise 

here  because  for  y  -  0,  the  normal  modes  of  the  system  decompose 

into  three  branches  with  the  modes  for  any  one  branch  involving 

displacements  of  the  particles  only  along  one  of  the  three 

2 

crystal  axes.  The  modes  associated  with  are  those  with 
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displacement  vectors  perpendicular  to  the  chains,  in  the  1 
direction.  If  y=  0,  the  forces  resisting  displacements  in 
the  1  direction  for  any  given  particle  arise  from  the  nearest 
neighbors  in  the  same  chain  and  from  the  nearest  neighbors 
in  those  adjacent  chains  which  lie  in  the  1,  3  plane  through 
this  particle.  The  essential  point  is  that  for  y  =  0  there 
is  no  coupling  between  the  displacements  of  this  particle  in 
the  1  direction  with  the  displacements  in  the  1  direction  by 
particles  in  adjacent  1,  3  planes.  Adjacent  1,  3  planes 
vibrate  independently  in  the  1  direction  and  so  for  this 
branch  with  y  =  0,  the  behavior  of  the  system  is  really  that 
of  a  two  dimensional  array  of  chains  rather  than  what  might 
be  described  as  a  typical  three  dimensional  array.  It  is 
this  feature  of  the  rectangular  lattice  which  is  responsible 
for  the  very  weak  dependence  of  on  q>2  mentioned  above. 

Similar  arguments  apply  to  the  branch  represented 

p 

by  v2,  except  that  the  role  of  the  1  and  2  directions  are  inter- 

2 

changed.  For  the  third  branch,  v^,  the  displacements  are  along 
the  chain.  The  dominant  feature  is  the  strong  central  force 
between  nearest  neighbors  in  the  chain*  The  nearest  neighbor 
interactions  between  chains  do  not  appear  in  because  the  cen¬ 
tral  forces  between  such  neighbors  act  only  perpendicular  to  the 
chain.  Any  coupling  between  chains  can  arise  only  from  the 
second  nearest  neighbors  and  therefore  is  very  weak.  That  the 
coupling  between  chains  is  relatively  weak  would  be  typical  of 
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most  chain  polymers  because  of  the  very  strong  central  forces  be¬ 
tween  adjacent  particles  in  the  chain  as  compared  with  the  forces 
between  chains.  However,  if  the  geometry  of  the  lattice  were 
such  that  nearest  neighbors  in  adjacent  chains  did  not  lie  in  a 
plane  perpendicular  to  the  direction  of  the  chain  then  the 
coupling  between  chains  would  be  proportional  to  the  nearest 
neighbor  force  constant  (a)  instead  of  the  second  nearest 
neighbor  force  constant  (y)» 

Our  purpose  here  is  to  analyze  specifically  the  model  as 
proposed  by  S-H  but  it  is  clear  that  the  peculiar  properties  of 
the  rectangular  lattice  are  closely  connected  with  the  choice 
of  y  «  a  and  that  if  y  were  not  chosen  small  compared  with  o, 
the  system  would  be  more  typical  of  chain  polymer  crystals  in 
general, 

2,  Vibration  Spectrum. 

The  first  problem  to  consider  is  that  of  analyzing  the 
vibration  spectrum  and  the  shape  of  the  surfaces  of  constant 
frequency  in  the  wave  number  space  <P2>  ^3^  From  S-H 

Eq,  (2,8)  we  note  that  and  Vg  differ  only  in  that  q^  and  q>2 
are  interchanged.  If  we  consider  the  entire  range  0  < 
q>2,  9^  <  it*  then  the  frequency  spectra  of  the  first  two  brandies 
are  identical.  The  third  branch  is  exactly  the  same  as  that  of 
a  cubic  lattice  which  has  already  been  considered  in  some  detail  • 
The  spectrum  of  the  third  branch  is  plotted  in  Fig,  5  for  the 
values  of  the  force  constants  used  by  S-H,  Only  the  spectrum 
of  needs  special  consideration. 
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The  frequencies  for  this  branch  are  given  by  S-H,  Eq,  2*8 


=  (a/2)(l  -  cos  9^)  +  y(2  -  cos  cp1  cos  cp2 

p 

-  cos  9^  cos  +  x(l  -  cos  9^)  (1) 

To  analyze  the  shape  of  the  surfaces  of  constant 
frequency  in  the  (9^,  92,  9  )  space,  we  first  inspect  the 
behavior  for  small  vj".  For  <  y  «  1,  it  is  necessary  that 
the  first  and  last  terms  of  Eq.  (1)  be  small,  of  order  y  °r  less 
with  y  «a,  x.  This,  in  turn,  requires  that  91  and  9^  both  be 
small  and  so  we  can  expand  in  powers  of  9^  and  9^  to  obtain 

~  [(ccA)  +  Cy/2)(1  +  cos  92)]92  +  Yd  “  cos  ?2) 

+  Cy/2)9^  +  (xA)9^  for  <q(Y)  (2) 


This  includes  all  terms  in  the  expansion  that  are  necessary  to 

p 

determine  the  qualitative  behavior  for  0  <  v-j_  <  y. 

For  sufficiently  small  v^,  namely  for  «  y2AH~10 
it  is  clear  that  92  must  also  be  small  and  that 
(y/2)92  »  (xA)9^  so  that  we  obtain  (S-H,  Eq.  3,4) 

vf  -  [(aA)  yD<pf  +  (y/2)?2  +  (Y/2)92  for  0  <  «^Ax 

(3) 

For  this  lowest  range  of  v^,  the  surfaces  of  constant  v 
are  flat  ellipses  with  axes  v^[(aA)  +  y]“^,  v^(y/2)"^  and 

v^(y/2)“‘2'  in  the  1,  2  and  3  directions  respectively.  The  first 
of  these  is  much  shorter  than  the  other  two  for  a  »  y,  see 
Fig,  1,  The  fact  that  the  surfaces  are  ellipoidal  for  sufficient¬ 
ly  small  v  is  a  feature  characteristic  of  many  three  dimensional 
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systems* 

As  v|  increases  past  y2Ax,  9 2  remains  small  compared 
with  one  hut  the  term  (xA)cp^  comes  into  prominence  and  quickly 
'overpowers  the  term  (y/2)<p^  because  of  the  assumption  x  >  >  y  • 
Whereas  the  <p^  and  <p2  intercepts  of  the  surfaces  of  constant  v 
continue  to  increase  nearly  proportional  to  v,  the  intercept 
increases  at  a  slower  rate  approaching  the  limit  (v^)z(xA)  • 

The  surfaces  become  elongated  in  the  cp^  direction  as  shown  in 
Fig.  1  and  we  see  that  (S-H  Eq.  3,6)  except  for  the  very  small 
values  of  9^, 


v|  -  [(oA)  +  Y^f  +  (Y/2)92  +  (XA)<?3  f*°r  Y2/lfH<<v1  «  Y 

A) 

This  form  of  reflects  the  influence  of  the  strong  valence 
forces  with  a  dependence  characterized  by  a  fourth  power 
term  rather  than  the  more  familiar  quadradic  dependence  as  in 


Eq,  3* 

p 

As  v£  approaches  the  value  Y»  "the  expansion  in  powers 
of  <p2  breaks  down  and  we  must  go  back  to  Eq.  (2)  which  tells 
us'  that  for  v2  =  2y,  the  surfaces  of  constant  v  pass  through  a 
corner  of  the  Brillouin  zone  at  (0,  it,  0)  creating  a  saddle 
point  such  as  discussed  by  van  Hove^*.  The  surface  v2  s  2y 
is  one  of  the  surfaces  shown  in  Fig.  2.  .  However,  we  shall  pass 
over  this  critical  range  of  v  for  the  moment. 

There  follows  another  region  of  frequencies 
Y  «  v£  «  a  in  which  the  y  terms  of  Eq,  (1)  are  of  only  minor 
importance*  As  can  be  seen  from  Fig,  2,  the  surfaces  of  con¬ 
stant  frequency  are  nearly  independent  of  q>2  (cylinders)  with 
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the  y  terms  only  causing  a  ripple  in  the  surface.  In  this 
range  of  v^,  still  small  and  except  for  the  very  smallest 

values  of  9^, 

v|  -  (a/2) (1  -  cos  9^)  +  (h/4)9^  fo r  y  «  v2  <  a 

(5) 

The  extent  to  which  such  a  range  of  frequencies  exists 
depends  upon  the  validity  of  the  condition  a  »  y.  Such  a 
condition  is  a  consequence  of  the  rectangular  geometry,  not 
simply  of  the  fibrous  structure. 

As  approaches  the  value  a,  the  surfaces  come  close 
to  an  edge  of  the  Brillouin  zone  and  the  detailed  shape  near 
this  edge  becomes  sensitive  to  the  y  terms  of  Eq.  1,  We  shall 
pass  over  this  critical  range  also  for  the  moment  and  consider 
still  larger  values  of 

For  a  «  v2  «  4*,  the  a  terms  of  Eq.  (1)  as  well  as 
the  y  terms  become  relatively  unimportant  but  9^  no  longer 
small.  v|  is  nearly  independent  of  both  9-^  and  92}  the  surfaces 
of  constant  are  nearly  planes  with 

v2  -  h (1  -  cos  9^)2  for  a  «  v2  <  4*  (6) 

As  approaches  the  value  4k,  the  surfaces  approach  a 
face  of  the  zone  and  the  dependence  of  on  a  and  y  once  again 
comes  into  prominence.  Figure  2  shows  some  of  these  surfaces 
of  constant  v,  particularly  those  which  pass  near  and  through 
critical  points  of  the  Brillouin  zone.  Because  of  the  very 
small  value  of  y,  some  of  the  de  °.ils  of  the  surfaces  can  not 
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be  shown  on  a  figure  of  this  scale  particularly  the  behavior  of 
the  surfaces  passing  near  the  two  corners  (it,  0)  and  (it,  0,  0). 

The  critical  values  of  the  frequency  which  were  passed 
over  in  the  above  description  give  rise  to  singularities  in  the 
frequency  distribution  to  be  discussed  in  the  next  section. 

These  are  associated  with  stationary  values  of  vf  considered  as 


a  function  of  q^,  <p2  and  cp^  and  can  be  found  analytically  from 

p 

Eq.  (1)  by  determining  those  values  of  (q>^,  <p2>  93)  and  vx  for 


which  dv|/a<p  =  0,  3  =  1,  2,  3.  These  critical  values  are  listed  • 


in  Table  I. 


From  the  signs  of  the  second  derivatives,  one  also 
determines  whether  any  such  point  represents  a  maximum,  minimum 
or  a  saddle  point  of  v|.  The  saddle  points  are  decomposed  into 
two  categories  denoted  by  S.P.I  and  S.P.II.  For  S,P,I,  the 
surfaces  of  constant  v  are  two  sheeted  for  frequencies  below 
the  critical  frequency  vc  and  one  sheeted  above  the  critical 
frequency,  S.P.II  differs  only  in  that  the  order  is  reversed; 
the  surfaces  are  two  sheeted  for  frequencies  above  the  critical 
frequency  and  one  sheeted  below.  This  distinction  shows  up 
in  the  frequency  distribution  g(v)  in  that  for  type  I  saddle 
points  there  is  a  verticle  tangent  singularity  as  v— ^  vc  from 
the  left  whereas  for  type  II,  the  verticle  tangent  singularity 

appears  as  v-^  v  from  the  right. 

0 
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Table  I*  The  critical  points  of  the  vibration  spectrum  for 
are  listed  below.  The  first  column  gives  the  point 
(cp^j  cp2»  ¥3)  wave  vector  space  and  column  two  the 

frequency  at  those  points.  The  third  column  gives  the  type 
of  singularity.  Max  denotes  local  maximum,  min,  minimum 
whereas  S.P.I  and  S.P.II  represent  two  types  of  saddle  points. 


■ 

92>  93) 

vf 

Type 

(0,  0,  0) 

0 

min. 

O 

O 

2r 

S.P.I 

(it,  x,  cos**^(l  -  y/2*)) 

a  +  2y  -  y2Ax 

S.P.I 

(x,  x,  0) 

a  +  2y 

S.P.II 

(x,  0,  cos^Cl  -  y/2h)) 

0  +  4y  -  Y2Ax 

S.P.II 

r~\ 

O 

O 

w 

a  +  4-  y 

max. 

/•% 

0 

0 

4*  +  2y 

S.P.I 

(0,  K,  lO 

4-x  +  4-y 

S.P.II 

(xj  Xj  h) 

4*  +  a 

S.P.II 

(x»  0,x) 

4x  +  a  +  2  y 

max. 
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It  is  interesting  to  note  that  there  are  more  than  the 
minimum  number  of  critical  points  (10  instead  of  8).  In  addi¬ 
tion  to  the  expected  singularities  at  corners  of  the  Brillouin 
zone,  there  are  also  two  saddle  points  at  (it,  0,  cos  (I-y/2*)) 
and  at  (it,  n,  cos“^(l  -  y/2k)).  Furthermore  the  singularity 
at  the  corner  (it,  0,  0)  is  a  local  maximum. 

Two  of  these  singularities  lie  very  close  to  the  corners 
with  <P2~(yA)^  ^  ^  The  frequencies  at  these  points  are 
extremely  close  to  the  corresponding  frequencies  at  the  corner, 
the  values  of  v2  differing  by  only  y^Ax  ~  10“^.  Such  details 
are  poorly  resolved  in  Fig,  2  and  do  not  appear  at  all  in  the 
spectrum  plotted  in  Fig.  4,  The  fact  that  there  are  two 
singularities  with  frequencies  very  close  together  is  of  little 
physical  importance  in  Itself;  it  is  of  interest  here  mainly 


because  we  shall  note  in  Sec,  4  that  as  the  two  singularities 
coalesce  into  one  for  y*~^°»  a  singularity  in  the  distribution 


of  a  slightly  different  type  results. 
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and 

0  <  9±  £  *  1  s  1,  2,  3 

It  is  possible,  however,  to  obtain  quantitatively  accurate 
approximations  to  g^Cv)  over  certain  ranges  of  v  based  upon 
the  approximations  used  in  Eqs,  (2)  and  (6). 

By  substituting  Eq.  (3)  into  Eq,  (7),  one  finds  very 
easily  that^* 

g^(v)  ~  2it“2Y~‘1'(a  +  v2<£v2  for  0  <  v2  «  y2/1™  (8) 

The  quadratic  dependence  of  g-^Cv)  on  v  is  typical  of  all  three 
dimensional  systems  and,  as  we  shall  see,  gives  the  usual 
Debye  T3  -law  for  the  specific  heat.  The  distinguishing  feature 
of  Eq,  (8)  is  the  fact  that  this  relation  holds  only  over  a 
very  short  range  of  v  as  compared  with  non-fibrous  crystals. 

From  Eq.  (^f)  we  find  that 

gx(v)  ~  2it”2Y’^(a  +  4Y)~V'*v3/2(£V3/2  for  y2A*  «  ^  «  Y 

(9) 

As  noted  also  in  Eq.  (4),  this  form  arises  because  at 
these  frequencies  the  valence  forces  have  become  the  dominant 
influence  in  the  <p^  direction.  The  v^2  dependence  is  a 
feature  that  is  expected  to  be  more  or  less  typical  of  chain 
polymer  crystals. 

Actually  by  combining  terms  in  Eq.  (3)  and  (4),  it  is 

possible  to  obtain  an  approximation  valid  over  the  entire  range 
o 

0  <  v  «  y  for  which  the  above  are  just  limiting  cases. 
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g-^v)  -  2^n“2(a  +  4y)”^y~^v  -^[(y/*)2  +  *+v2/x]^  -  y/x}^ 

for  0  <  v2  «  y  (10) 

Equation  (10)  is  plotted  as  the  solid  curve  of  Fig*  3  and  the 

two  limiting  equations,  Eqs,  (8),  (9)  are  shown  by  the  broken 

lines.  The  range  of  frequencies  involved  here  is  very  low  and 

the  whole  of  Fig,  3  is  a  greatly  magnified  (about  40  to  1)  view 

of  the  small  rectangle  near  the  origin  of  Fig.  4-. 

2 

If  we  pass  over  the.  region  around  v  =  2y  for  the  moment 
and  go  to  slightly  higher  frequencies,  we  come  to  a  range  where 
Eq.  (5)  holds.  Despite  its  simple  form,  Eq.  (5)  does  not  yield 
elementary  integrals,  but  leads  to  the  form 

gx(v)  •  u"2x^^2  dp-Jv2  -  a  sin2(91/2)]“3^  (11a) 

4) 

2  2 

in  which  9  is  the  value  of  9.  for  which  v  =  a  sin  (9,/2). 
mi  -*• 

If  we  let  y  =  av"2sin2(9^/2),  Eq.  (11a)  gives 

1 

g_(v)  ~2“^7f2,f*(a  +  dy  y~k(l  -  y)"^(l  -  v2y/a)~^ 

1  JO 

(lib) 

and  if  we  restrict  v  so  that  v  «  a,  the  integral  can  be 
evaluated  in  a  power  series  in  v2/a  of  which  the  leading  term  is 

gx(v)  -  |  iT5/W*[r(ij)]2  »V  for  2y  «  V2«  o  (12) 

T  is  the  gamma  function  which  can  also  be  expressed  in  terms 
of  complete  elliptic  integrals  through  the  relation 

cn^)]2  =  4  vests'*). 

We  could  also  have  obtained  Eq,  (12),  by  restricting  the 
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2 

range  of  9^  in  Eq.  (5)  and  replacing  cos  <p^  by  1  -  9^2* 

The  v"^2  dependence  results  from  having  cylindrical  wave 

1  /2 

surfaces  with  one  intercept  proportional  to  v  '  and  the  other 
proportional  to  v  and  as  noted  in  Eq*  (5),  this  feature  is 
mainly  due  to  the  rectangular  geometry.  With  a  value  of  a/2y 
of  only  about  10,  the  existence  of  a  well  defined  frequency 

p 

range  with  2y  «  v  «  a  may  be  questioned  even  for  the  values 
of  the  constants  used  here,  If,  however,  we  could  make  this 
ratio  appreciably  larger  such  a  range  would  certainly  exist. 

A  plot  of  Eq.  (12)  is  shown  in  Fig.  4  as  one  of  the  broken 
lines.  It  is  seen  to  give  only  a  fair  approximation  over  a 
quite  limited  range  and  for  these  values  of  the  constants,  the 
more  accurate  curve  shown  by  the  solid  line  would  hardly  be 
recognized  as  being  proportional  to  v1/2. 

If  we  consider  next  the  region  of  still  higher  frequencies 
where  a  «  v|  <  4k,  we  can  use  Eq.  (6)  which  again  yields  a 
simple  expression  for  g^(v) 

g;L(v)  -  2~*,f1x"*[l  -  (v2/4H)*r*v~^  for  -  *a  «  v2  <  4k  (13) 

which  actually  gives  a  fairly  good  approximation  over  most  of 
the  range  from  a  to  4k*  This  equation  is  not  shown  in  Fig*  4 
but  the  solid  curve  shown  is  obtained  from  a  power  series 
expansion  in(a/K)  that  converges  quite  rapidly  over  most  of 
this  range.  Equation  (13)  is  the  lowest  term  of  this  expansion. 

It  is  interesting  to  note  that  from  Eq.  (13),  we  can 
obtain  another  limiting  form 
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g  (v)  -  2“V1x"V^  «  v"^  for  <£  «  v2  «  4x  (14) 

1 

Again  for  the  choice  of  constants  used  here  with  4*/a  *  10, 
the  existence  of  any  such  frequency  range  is  questionable*  A 
plot  of  Eq,  (14)  is  shown  in  Fig.  4  where  one  may  see  that  this 
approximation  is  not  particularly  good.  The  v"^  dependence 
of  g(v)  does,  however,  represent  a  limiting  situation  in  which 
4# /a  >  >1  and  would  indeed  be  recognized  if  we  could  make  this 
ratio  appreciably  larger  than  just  10. 

This  frequency  dependence  is  not  simply  a  result  of 
geometry*  Similar  results  would  also  be  obtained  for  x  »  a,  y 
without  any  assumptions  regarding  the  relative  size  of  a  and  y  • 
It  is  a  consequence  of  the  choice  of  the  valence  forces  as  the 
dominant  influence  for  the  modes  vibrating  perpendicular  to  the 
chain  and  is  likely  to  be  typical  of  any  chain  polymer  crystal. 

It  is  worth  noting  perhaps  at  this  time  that  some  of  the 
approximations  discussed  so  far  work  very  well  and  others  are 
rather  crude  even  though  the  assumption  a/2y  ~  10  »  1  and 
4*/ a  ~  10  »1  seems  quite  reasonable.  The  reason  for  this  is 
that  the  errors  are  in  some  cases  proportional  to  (2Y/<x)^  or 
(cc/4h)£  ~  ^  or  even  (a/4x)^~  rather  than  (a/4x)  •*'  For 
example,  we  may  say  that  is  small  in  Eq,  (5)  when 
v2  ~  a  <<  x  but  actually  Is  order  (ba/*)*  which  really 
is  not  very  small  $  in  fact  it  is  very  nearly  1*  It  is  not 
surprising  then  that  the  approximations  in  Fig*  4  are  not 
accurate  to  the  order  of  10#  as  might  be  expected  from  a  hasty 
judgment.  Indeed  such  things  as  effective  Debye  temperatures 
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and  other  qualitative  measures  of  the  sizes  of  physical  con¬ 
stants  are  really  proportional  to  frequencies  which  in  turn 
are  proportional  to  square  roots  of  force  constants  rather  than 
the  force  constants  themselves. 

On  the  basis  of  the  analysis  so  far,  we  see,  however, 
that  for  an  extreme  situation  in  which  y  «  a  «  k,  the  frequency 
distribution  g-^Cv )  varies  as  v2  at  the  lowest  frequencies, 
gradually  changes  over  to  a  v^/2  dependence  starting  at  v2  in 
the  vicinity  of  y^/4x,  later  goes  over  to  a  v-^2  dependence 
when  v  becomes  larger  than  y  and  finally  varies  as  v"  after 
v2  becomes  larger  than  a  but  still  small  compared  with  4x*  The 
effects  of  this  on  the  specific  heat  will  be  shown  in  Sec*  5* 

The  problem  of  actually  constructing  an  accurate  graph 
of  g^(v)  for  the  entire  range  of  v,  is  a  rather  tedious  job* 

For  each  of  the  frequency  ranges  considered  above,  the  approxi¬ 
mate  formulae  for  g-^(v)  can  be  used  as  first  approximations  and 
more  accurate  expressions  obtained  by  expansions  in  powers  of 
y/a,  a/x  (fractional  powers  in  most  cases)  or  whatever  may  be 
appropriate,  depending  upon  what  has  been  neglected  in  Eq*  (1), 
to  arrive  at  the  lowest  approximation.  Such  a  procedure  was 
used  to  obtain  the  solid  curve  of  Fig.  b  in  the  appropriate 
ranges.  The  formulae  and  detailed  procedures  for  this  are 
f a Jjtfly  straightforward  and  not  sufficiently  interesting  to 
warrant  presentation  here.  Suffice  it  to  note  that  this 
procedure  is  not  a  satisfactory  one  for  obtaining  g^(v)  near 
its  singularities  where  the  expansions  will  converge  very  slowly 
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If  at  all.  The  regions  considered  by  this  method  are  non¬ 
overlapping  and  give  little  information  regarding  g-j_(v)  for  v2 
near  the  values  2 y,  a  or  4x, 

Singular  points 

Van  Hove^*  has  described  the  analytic  form  of  the  singular 
points  for  any  quite  general  two  or  three  dimensional  crystalline 
lattice.'  This  description  can  indeed  be  applied  to  the  present 
problem  but  it  really  does  not  tell  the  whole  story.  Highly 
anisotropic  crystals,  in  general,  and  the  present  model  in 
particular  with  strong  valence  forces  give  rise  to  anomalous 
behaviors  near  many  of  the  singular  points  as  illustrated,  for 
example,  in  Pig.  4  at  the  point  y2  *=  o»  This  singularity, 
though  actually  having  a  verticle  tangent  only  on  one  side  of 
the  singular  point,  appears  more  like  a  true  cusp  with  a  verticle 
tangent  on  both  sides  of  the  singularity. 

We  have  already  seen  that  the  strong  valence  forces 
cause  unusual  behavior  near  the  origin,  the  minimum  point  of 
the  spectrum,  which  is  in  a  certain  sense  one  of  the  singu¬ 
larities.  Also  we  have  seen  that  for  v2  near  a  +  2y  and  also 
for  v2  near  a  +  4-y,  there  are  two  singularities  so  close 
together  that  it  would  be  impossible  to  resolve  them  on  a  scale 
of  frequencies  such  as  that  used  in  Pig. 

Table  I  lists  all  the  singular  points.  By  expansion  of 
v  as  given  by  Eq.  (1)  about  the  critical  values  of  9^,  92f  9^ 
one  can  easily  deduce  in  accordance  with  van  Hove  that  near  any 
critical  frequency  vq,  except  v  »  0,  g(v)  has  the  general  form 
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g(  v)  =  a  +  b  Real  [  +  (  v  -  v)]^+0(v  -  v)  (15) 

c  c 

The  +  sign  inside  the  square  root  and  the  sign  of  b  are 
chosen  in  accordance  with  the  four  types  of  singularities  as 
listed  in  Table  I.  b  is  negative  for  S.P.I  and  S.P.II,  positive 
for  max*  and  min*  The  +  sign  is  chosen  for  S.P.II  and  min., 
the  -  sign  for  S.P.I  and  max.  The  values  of  the  constant  b  is 
determined  easily  from  values  of  the  second  derivative  of  v2 
with  respect  to  the  <p's  at  the  critical  <p*s.  The  value  of 
the  constant  a  does  not  usually  depend  only  upon  the  local 

p 

properties  of  v  near  the  critical  9*3,  however,  and  is  not  so 
easily  determined  explicitly. 

The  origin  v  =  o  does  not  give  a  singularity  in  g(v) 
but  there  is  a  vertical  tangent  singularity  in 

G(v2)  =  (2V)"1  g(v)  (16) 

the  distribution  of  values  of  v2.  In  place  of  Eq.  (15)  we 
find  for  G(v2)  the  analogous  form 

G(v2)  =  A  +  B  Real  [  +  (v2  -  v2)]*  +0(v2  -  v2)  (1?) 

at  every  critical  frequency  vQ  including  v  =  0, 

The  behavior  of  G(v2)  near  the  origin  indicates  to  a 
oertain  extent  the  type  and  cause  of  anomalies  in  G(v2)  and 
g(v)  at  some  of  the  other  critical  points, 

p 

From  Eq.  (8)  and  Eq.  (16),  we  see  that  G(v  ),  considered 

2 

as  a  function  of  v  ,  starts  out  at  v  =  0  with 
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G<v2)£v  =  (v2)^  for  0  <  v2  «  y^Ax  (18) 

in  agreement  with  the  general  form  Eq.  (17)  with  A  =  0  (as 

it  must  be  for  an  absolute  minimum). 

v2  has  a  complete  range  of  order  1  whereas  Eq.  (18)  holds 
only  for  v2  «  ~  10“^,  This  behavior  is  confined  to 

such  a  small  range  that  it  would  not  be  resolved  on  any  graph 
showing  the  whole  spectrum  G(v2)  vs,  v2;  it  did  not  even  show 
on  the  graph  of  g(v)  vs.  v  which  is  relatively  much  more 
sensitive  to  the  behavior  near  v  =  0. 

As  soon  as  the  valence  forces  come  into  prominence  with 
the  characteristic  dependence  of  v2  on  9^,  g(v)<£v^2  and 

G(v2)<£v^  «  (v2)^  for  y^A*  «  v2  «  2y  (19) 

The  singularity  has  changed  from  a  square  root  type  to  a  fourth 

root.  Although  this  formula  is  also  confined  to  a  very  small 
scale  of  v2  values,  v2  «  k  x  10”^,  it  is  this  feature  which 
is  of  particular  interest  for  it  is  one  of  several  possible 
types  of  singularities  that  arise  when  v  has  a  stationary  value 
with  a  local  dependence  upon  the  wave  numbers  that  is  quadratic 
in  some  directions  (9^  and  92)  but  primarily  fourth  power  in 
other  directions  (9^).  Ibis  feature  is  characteristic  of  all 
the  singularities  that  lie  on  or  near  the  face  9^  s  0  in  this 
particular  problem.  It  is  a  feature  that  is  not  peculiar  to 
this  particular  model  but  will  arise  in  most  any  model  having 
strong  valence  forces  with  the  characteristic  fourth  power 
dependence  of  v2  on  one  or  more  of  the  wave  numbers. 
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The  fourth  root  singularity  is  typical  of  any  system 
2 

in  which  v  behaves  locally  as 

y2  s  vc  ^  a(<p!  "  9lc)2+  b(<p2  "  ?2c)2  +  c(<?3  "  ^30^  (20) 

in  which  vc  is  a  critical  frequency,  a,  b,  c  are  positive 
constants  and  (<Plc>  <P2c»  ^3c^  a  cri^ca^  wave  vector. 

Because  of  the  rectangular  geometry  and  the  fact  that 
this  model  degenerates  into  a  two  dimensional  problem  with 
Y — ^0,  we  have  an  illustration  here  of  still  another  type  of 
singularity.  Equation  (5)  for  ^  «  1  represents  a  two  dimen- 
sioanl  form  of  Eq.  (20),  b  =  0,  and  as  we  have  already  seen 
g(v>$v^2  and 

G(v2)<fcv~k  s  (v2)-±  for  2y  <<  v2  «  a  (21) 

The  negative  fourth  root  singularity  is  typical  of  the  two 
dimensional  form  of  Eq.  (20). 

On  still  a  larger  scale  of  v2,  Eq.  (6)  with  93  «  1 
represents  a  one  dimensional  form  of  Eq.  (20)  with  a  =  b  =  0 
and  g(v)<,v^ 

G(v2)gy"3^2  =  (v2r3A  for  a  «  v2  «  4x  (22) 

The  negative  three-fourths  powers  singularity  is  representative 

of  the  one  dimensional  form  of  a  minimum  having  a  fourth  power 

dependence  of  vc  on  <p_. 

J  2 

If  now  we  go  on  to  the  singularity  at  v  «  2y,  we  do 
indeed  see  that  in  the  immediate  vicinity  of  ( 09  n,  0),  there  is 
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a  singularity  of  standard  type  (S.P.I)  but  since  it  is  on  the 
face  s  the  dependence  of  v2  on  9^  is  to  a  large  extent 
dominated  by  the  fourth  power  term*  v2  behaves  approximately 
as  given  in  Eq,  (20)  but  with  a  negative  value  of  b.  We  have 
a  "singular"  type  of  saddle  point. 

On  the  scale  with  which  Fig.  b  is  drawn,  only  the  behavior 
typical  of  the  singular  type  of  saddle  point  can  be  shown.  We 
have  included  in  the  appendix  a  discussion  of  all  the  possible 
types  of  singularities  that  might  arise  from  valence  forces. 

It  is  shown  there  that  this  type  of  saddle  point  gives  rise  to 
a  fourth  root  verticle  tangent  singularity  in  g(v)  on  both 
sides  of  vQ  (not  necessarily  symmetric  about  vc,  however). 

This  explains  why  g(v)  appears  to  have  a  true  cusp  at  v  =  2y, 

It  has  been  possible  to  actually  obtain,  explicitly, 
accurate  analytic  approximations  for  g(v)  everywhere  including 
all  the  singular  points.  Unfortunately  the  derivation  of  these 
expressions  involves  using  slightly  different  techniques  at 
each  of  the  singular  points  depending  upon  which  of  the  many 
variables  involved  in  the  integrations  may  be  small  for  any 
particular  value  of  v.  The  results  are  not  exactly  trivial  and 
in  most  cases  involve  combinations  of  T-functions.  It  would 
not  be  possible  to  describe  all  the  various  approximations  here. 
Mainly  as  an  illustration  of  the  form  of  a  typical  result  we 
give  below  an  approximation  for  g(v)  which  is  fairly  accurate 
in  a  small  region  near  the  point  v2  «  2y 
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g(v)  =  fl  -  Uf,A,  Heal  f-it 

2(a  +  2Y)*k4(2x)H  [r(lA)  t  l  2Y 
for  10“3  -  y/16  x  «j-  1  +  v/  y/2y  I  «  1 

This  formula  does  not  apply  at  the  singular  point 
itself  because,  strictly  speaking,  there  is  no  verticle  tangent 
on  both  sidee  of  the  singularity  but  it  would  be  necessary  to 
draw  a  graph  that  shows  frequency  differences  of  order 
v/2y  Y/16  h  ~  10"4  to  see  this. 

The  behavior  near  v2  ~  a  is  even  more  unusual.  Near 

p 

the  point  (it,  it,  0),  v  -  a  +  2yf  there  are  two  singularities, 

p 

This  results  because  near  this  point  the  dependence  of  v  on 
<p^  is  represented  by  a  negative  but  small  (proportional  to  y) 
quadratic  term  plus  a  positive  but  large  (proportional  to  x) 
fourth  power  term.  As  a  function  of  9^,  v2  ^as  a  minimum  result 
ing  from  the  competing  influences  of  the  quadratic  and  fourth 

p 

power  terms.  In  addition  to  the  stationary  value  of  v  at 
9^  =  0  there  is  a  second  stationary  value  a  short  distance  away 
from  the  corner.  This  anomaly  is  at  least  indicated  in  Fig.  2 
but  the  double  singularity  cannot  be  resolved  in  Fig,  Figure 
V  shows  only  the  influence  of  the  fourth  power  9^  terms  but  not 
the  quadratic  terms.  Disregarding  this  fine  structure,  the 
behavior  of  v2  is  described  by  an  equation  of  the  type  Eq,  (20) 
with  a  <  0  and  b,  c  >  0.  This  gives  another  singular  saddle 
point  of  exactly  the  same  type  as  at  v2  =  2y  with  approximately 
verticle  tangents  for  g(v)  on  both  sides  of  the  point 
v2  e  o  +  2f. 


f) 
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At  the  point  (it,  0,  0),  v2  ~  a  +  4y,  we  have  a  similar 
situation  with  a  pair  of  singular  points  due  to  the  same  sort 
of  dependence  of  v2  on  9^.  The  exact  picture  is  somewhat  more 
complicated  than  at  the  point  (71,  7t,  0)  because  of  the  creation 
at  the  corner  of  a  local  maximum  and  an  "island"  by  the  saddle 
point  a  short  distance  away  from  the  corner.  Figure  4  is,  how¬ 
ever,  again  drawn  on  too  coarse  a  scale  to  show  this  and  near 
this  corner  the  frequencies  are  given  approximately  by  Eq.  (20) 

with  a,  b  <  0  and  c  >  0.  Although  this  is  also  a  singular 

2 

saddle  point,  it  is  not  of  the  same  type  as  those  at  v  =  2y 

and  v2  =  a  +  2y  by  virtue  of  the  fact  that  both  quadratic  terms 

(a  and  b)  are  negative  at  v2  =  a  +  whereas  at  the  other  two 

frequencies  only  one  of  these  terms  is  negative.  It  is  shown 

in  the  appendix  that  this  new  type  of  saddle  point  gives  rise 

to  an  infinite  discontinuity  in  slope  with  a  finite  slope  on 

one  side  (in  this  case  the  low  frequency  side  of  v  )  and  a 

c 

vertical  tangent  on  the  other  side.  In  this  respect  it  is  the 
same  sort  of  behavior  as  obtained  from  a  saddle  point  of  stand¬ 
ard  type  but  the  verticle  tangent  arises  in  the  present  instance 
from  a  fourth  root  rather  than  a  square  root  singularity. 

The  fact  that  for  y  <  <  1  all  singular  points  at 
v2  =  a  +  2y  and  v2  =  a  +  *ty  lie  very  close  together  indicates 
that  on  a  coarse  scale  of  frequencies  in  which  frequency  differ¬ 
ences  of  order  y  are  not  resolved,  the  g(v)  curve  has  an 
approximate  singularity  of  still  another  kind  in  the  vicinity 
of  v2  =  a  •  This  is,  of  course,  due  to  the  fact  that  for 
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Y  — £  0,  the  problem  reduces  to  a  two  dimensional  problem  and 
as  a  two  dimensional  problem,  the  surfaces  of  constant  v  are 
all  cylinders  independent  of  g>2  the  line  (it,  q>2,  0)  a 

line  of  singularities# 

The  form  of  the  surfaces  of  constant  frequency  is 
represented  approximately  by  Sq.  (20)  with  b  =  0,  a  <  0,  c  >  0, 
a  singular  form  of  a  two  dimensional  saddle  point#  In  the 
appendix  it  is  shown  that  this  gives  rise  to  an  infinite  singu¬ 
larity  of  the  type 


g(v)  oC  Real  (v 

c 


2-*Cv  - 


for  vc  >  v 
for  vc  <  v 


Although  g^(v)  has  a  minus  one  fourth  power  singularity 
on  both  sides  of  vc,  g1(v)  Is  not  symmetric  about  vc#  Figure 
4  clearly  shows  the  tendency  to  approximate  this  two  dimensional 
behavior  for  v2  ~  a  • 

The  remaining  singularities  of  the  spectrum  all  occur 
on  the  face  a  x  and  are  of  a  more  conventional  variety 
although  we  can  again  see  the  effects  of  having  a  spectrum  that 
on  a  very  coarse  scale  tries  to  approximate  a  one  dimensional 
spectrum,  on  a  finer  scale  a  two  dimensional  spectrum  and 
finally  the  three  dimensional  spectrum  on  a  still  finer  scale# 

We  no  longer  have  a  local  fourth  power  dependence  of  v2  on  the 
wave  numbers  near  singularities  as  is  true  on  the  face  9^  55  0* 
The  curve  clearly  shows  the  usual  varieties  of  verticle 
tangent  singularities  characteristic  of  any  three  dimensional 
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O 

system  at  the  four  points  v  =4-x  +  2y,  4x  +  4y»  4-x  +  a  and 

4-x  +  a  +  2y  which  arise  from  the  stationary  values  of  v2  at 

the  four  corners  of  the  cube  on  the  faces  9^  =  On  a  scale 

in  which  frequency  differences  of  order  y  are  not  resolved,  the 

singularities  coalesce  in  pairs  at  the  approximate  points 

v2  =  4*  and  v2  =  4-x  +  a*  In  this  approximation  the  spectrum 

is  two  dimensional  with  the  usual  two  dimensional  singularities, 

a  logarithmic  singularity  at  v2  =  4-x  coming  from  the  edge 

(0,  9  ,  it)  of  Fig,  2  and  a  discontinuity  at  the  maximum  frequency 

4x  +  a  coming  from  the  edge  (it,  92,  it)  of  Fig,  2,  This  tendency 

of  Fig,  4-  to  approximate  this  behavior  is  clearly  indicated. 

On  a  still  coarser  scale  of  frequency,  all  four 
singularities  occur  at  approximately  the  same  frequency, 
v2  -  4x  and  the  spectrum  approximates  that  of  a  one  dimensional 
system  with  the  characteristic  recriprocal  square  root  singular¬ 
ity  in  g(v)  due  to  the  maximum  on. the  face  of  the  cube 
(9^,  92,  it) ,  This  tendency  is  illustrated  in  Fig,  4-  by  the  fact 
that  g^(v)  rises  quite  rapidly  as  v2— ^  x  on  a  coarse  scale  of 
frequency. 

The  fact  that  it  is  possible  actually  to  construct  a 
fairly  accurate  curve  for  g^(v)  throughout  the  whole  range  of 
v  is  to  a  large  extent  a  consequence  of  the  fact  that  the 
spectrum  is  really  a  distortion  of  a  one  dimensional  spectrum 
into  a  two  dimensional  spectrum  and  then  a  distortion  of  the 
two  dimensional  spectrum  into  a  three  dimensional  spectrum# 

No  attempt  has  been  made  to  analyze  carefully  the  size 
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of  the  errors  made  in  the  calculation  of  g^(v)  but  the  errors 
of  Pig,  4-  exceed  a  few  percent  only  at  a  few  places  at  most. 

The  curve  is  at  least  as  accurate  as  we  believed  anyone  would 
ever  care  to  know  it, 

5*  The  Complete  Spectrum 

As  noted  in  Sec.  2,  the  total  spectrum  consists  of  three 
branches  of  which  the  first  two  are  identical  and  have  been 
described  in  the  previous  sections.  The  third  branch  is  the 
same  as  that  of  a  simple  cubic  lattice  and  Fig,  5  is  drawn 
using  the  formulae  given  in  Ref.  4.  We  wish  merely  to  point 
out  a  few  features  of  this  spectrum  that  were  not  noted  there. 
The  spectrum  is  approximately  one  dimensional  and  in 
contrast  with  the  somewhat  artificial  application  intended  in 
Ref,  4,  this  appears  to  be  a  fairly  realistic  feature  of  the 
fibrous  crystals  although  this  is  perhaps  exaggerated  by  making 
y/p  ~  2  x  10~3 # 

g^(v)  starts  at  v  =  0  with  the  characteristic  v2 
dependence 

g3(v)  -  27f2Y~1(p  +  4Yr*v2Jlv2  for  v2  «  2Y  (23) 

There  are  singularities  at  v2  *=  2Y  and  4y  but  as  v2  becomes 
larger  g(v)  becomes  almost  equal  to  the  one  dimensional  spectrum 
for  which  g^ ( w)  is  a  constant  at  low  frequencies 

g3(v)  -  -  v°  for  <  v2  «  1  (24) 

These  two  limiting  forms  are  shown  in  Fig,  5  by  the  broken  line. 
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The  fact  that  y/j3  is  so  small  makes  these  quantitatively 
accurate  in  contrast  with  some  of  the  less  accurate  limiting 
forms  for  g^v). 

At  the  high  frequency  end  of  the  spectrum  ( v)  tries 
to  become  infinite  as  its  approximate  one  dimensional  form 
requires.  The  three  dimensional  form  requires  that  g(v)  remain 
finite  as  indeed  the  magnified  view  of  this  end  of  the  spectrum 
shows  in  Fig.  5* 

The  complete  spectrum  of  the  system  is  the  sum  of  the 
three  branches 

g(v)  =  ^g^v)  +  g^(  v)  ]/3 

Fig,  6  shows  the  resulting  distribution  by  the  solid  curve. 

The  broken  line  step  function  is  the  numerical  graph  obtained 
by  S-H. 

S-H  recognized  that  the  numerical  calculation  with  a 
rectangular  array  of  mesh  points  in  the  wave  number  space 
converged  very  slowly  for  computing  g^ (v).  This  is  true  not 
simply  because  of  the  odd  shape  of  g^(v)  but  because  of  the 
fact  that  is  nearly  independent  of  9^  and  9 2  and  the  mesh 
points  run  almost  parallel  with  the  surfaces  of  constant  v. 

Very  little  increase  in  accuracy  results  from  taking  a  mesh 
with  more  than  one  value  of  (9^,  92),  Tb©  accuracy  depends 
mainly  upon  the  number  of  values  of  93  that  are  used,  S-H  did 
use  a  finer  mesh  for  evaluating  g^( v)  than  for  g^(v)  but  the 
same  difficulties  also  arise  to  a  lesser  extent  in  computing 
the  central  portion  of  g^(v)  for  which  the  surfaces  of  constant 
v  are  also  nearly  planes. 
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This  probably  explains  the  apparent  systematic  error 
between  v  -  0.2  and  v  ~  0.6.  (The  solid  curve  is  accurate 
to  within  a  fraction  of  a  percent  in  most  of  this  region.)  At 
lower  frequencies  the  surfaces  of  constant  v  are  much  less  flat 
and  the  statistical  errors  in  the  numerical  scheme  should  be 
much  less  for  0  <  v  <  0,2,  The  agreement  is  indeed  very  good 
in  this  range. 


6,  The  Specific  Heat 

The  specific  heat  at  constant  volume  is  evaluated 
from  the  equation 

1 

r* 

G  =  3R 


VT 


dvg(v)f(vTm/T)  =  3R(T/Tm) 


dx  g(xT/Tm)f(x)  (25) 


in  which  v  is  still  the  relative  frequency  as  used  in  the 
previous  sections,  f(x)  is  the  Einstein  function 


and 


f (x)  =  x2ex(cx  -  l)’2 


T_  =  hv  /k,  x  =  vT  /T 
mm*  m 


(26) 


Tm  is  the  temperature  (analogous  to  the  Debye  temperature) 
associated  with  the  maximum  frequency  vm>  for  which  v  =  1, 

For  low  temperatures,  T  «  Tffl,  only  the  low  frequency 
modes,  v  «  1,  contribute  appreciably  to  C.  Since  the  function 
f(x)  decreases  rapidly  for  x  »  1  the  main  contribution  to 


Eq#  (25)  comes  from  values  of  x  of  order  1  and 

dx  g(xT/Tm)f (x)  for  T  «  Tm 


00 

A 


C  -  3R(T/Tm) 


(25a) 


0 


f 


%ntM 


a 


Nonr  562-08 


29 


If  for  small  values  of  v  =  xT/Tm,  g(v)  should  be  pro- 

^  ft 

portional  to  v  for  some  value  of  l,  i.e,,  g(v)  *=  A then 


C  -  3RA(T/Tm) 


l+l 


oo 

r* 


JL  ixi 

dx  x  f (x)  =  3HAr(>t  +  3)C(-t  +  2) (T/Tm) 

(27) 


oC  T 


t+1 


CC-t  +  2)  is  the  Riemann  zeta-f unction. 

It  is  convenient  to  consider  separately  the  contribution 
to  C  from  the  branches  g.^(v)  =  g2(v)  and  g^(v)  ^y  writinS 

C  =  (2CX  +  C3)/3  (28) 

for  which  is  the  contribution  to  C  from  the  branch  g1(v)  and 
the  contribution  from  g3(v). 

For  v2  «  Y2A*t  g^(v)  is  proportional  to  v2,  Eq.  (8), 
and  for  v2  «  y»  g^(v)  is  also  proportional  to  v2,  Eq.  (23), 

In  the  range  x  of  order  one,  v  =  xT/Tm  «  y/  if 

T  «  T  y/  v/^+k  and  so  we  find  that 
4Ric2(T/Tm)3 


C,  - 


1  5y(a  +  ky)t 


and  similarly 


t£.T^  for  T  <<T_  y/  y^fn  «  T 
i  m  rn 


(29) 


C.  ~ 


4Rn2(T/Tm)3 


<£  T3  for  T  «  Tm  Vy  «  Tj 


m 


(30) 


3  5Y(P  +  *+y)* 

C  as  well  as  its  two  components  individually  is 
proportional  to  for  sufficiently  low  T  in  agreement  with  the 
familiar  Debye  T^law  except  that,  in  the  present  casef  the 


9 
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behavior  extends  over  only  a  relatively  very  short  range  of  T; 
for  a  typical  three  dimensional  system  the  T^  behavior  is  limited 
only  by  T  «  Tm  i.e.,  for  temperatures  small  compared  with  the 
Debye  temperature. 

From  Eq,  (9)  we  see  that  g^(v)  is  proportional  to  v^2 
for  y2Ax  «  v2  «  y  •  Y/*  is  sufficiently  small,  the  v^/2 

p 

dependence  holds  over  almost  the  entire  range  v  «  y»  except 
for  a  very  small  region  near  v  =  0,  represents  an  average 
of  contributions  from  all  values  of  v  =  xT/Tm  with  x  =  0(1), 

If  TmY/\An  «  T  «  Tm\/Y  >  then  g(v)oCv^2  for  most  values 
of  v  that  contribute  to  C-^  and  we  find  that 


315H(T/Tm)5/\(7/2)  5/2 

8x3/2Y1/2(a  +  i¥y^l/2)ll/k 


for  TmY/  AT  «T«Tm  yft 

(3D 


5/2 

The  T  dependence  of  the  specific  heat  is  not  a  con- 

■ i/p 

sequence  of  geometry  and  the  dependence  was  one  of  the 

better  approximations  in  the  analysis  if  g^Cv).  At  least  a 
trend  in  this  direction  should  appear  for  any  fibrous  crystal. 
If  one  should  plot  vs.  T  on  a  log-log  scale,  the  curve 
should  have  a  slope  3  at  the  lowest  values  of  T  and  show  a 
range  of  T  over  which  the  curve  is  rather  flat  with  a  slope 
of  5/2, 


One  does  not,  of  course,  measure  experimentally  but 

only  C  and  in  this  range  of  T,  we  expect  to  be  proportional 

to  T^  still.  The  total  specific  heat  would  be  a  linear  com¬ 
bination  of  T^  and  T^/2  and  the  log-log  plot  of  C  would  have  a 
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slope  intermediate  between  3  and  5/2,  nearer  the  5/2,  however, 
because  at  these  temperatures  would  be  appreciably  larger 
than  Cy 

For  this  particular  model,  we  would  find  also,  by 
quite  similar  arguments,  a  range  of  T  for  which 


45R(T/Tm)3/2[r(lA)]2i(5/2)  V2  r-  • 

Cl - —  TTTT  Va -  *  T  for  ^  «  T<^“  Tm 

16,1  °  “  (32) 

1/2 

due  to  the  v  '  dependence  of  g^(v)  given  in  Eq.  (12),  This 
behavior  of  g^(v)  was  a  limiting  form  that  did  not  show  too 
well  for  the  force  constants  used  here  and  furthermore  it  was 


primarily  a  consequence  of  the  rectangular  geometry.  The 
specific  heat  being  an  average  of  contributions  from  a  wide 


range  of  v,  would  not  be  expected  to  show  this  behavior  even 

n  3/2 

as  accurately  as  g^(v)  showed  the  dependence.  The  T 

represents  only  a  trend  peculiar  to  the  geometry. 

Also  from  Eq.  (1^-)  which  gives  g^(v)  cC  we  obtain 

9R(T/Tm)V\(3/2)  1/2 

C1  '  2^2Bl/2KlA  *  T  for  T,  «  1  «  **  <33) 


This  is  expected  to  be  only  a  rather  crude  approximation  as 
was  Eq,  (Ik)  because  it  requires  sAk/cx  »  1,  This  dependence 


was  not  mainly  a  consequence  of  geometry,  however  and  we  can 
expect  this  sort  of  temperature  dependence  to  be  approximately 
true  for  any  fibrous  structure. 

Returning  to  the  behavior  of  we  see  also  from  Eq# 
(2k)  that 
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C3  ~  2  R*3<T/V  *  T  for  2  fy”  Tm  «  T  «  Tm  (3*0 

This  behavior  is  not  mainly  a  consequence  of  geometry  but  of 
the  fact  that  the  central  forces  in  the  chain  are  very  strong* 
In  a  crude  sort  of  way,  we  can  eliminate  the  effects 


of  the  geometry  by  making  y  of  the  same  order  of  magnitude  as 
a  •  The  features  to  be  expected  of  any  fibrous  structure  are 
that 


CrflC  T3 

for  T  <  <Tm  a/  v^+x 

<£  T5/2 

for  Tffi  a/  «  Tm  «  \/a  Tm 

Cx  oC  T1/2 

for  T  \/a  «  T  «  T 
m  mm 

and 

for  T  «  y/a  Im 

C3<C  T-3 

C3cC  T 

for  /a  T  «  I  «  T 
'  m  m 

In  which  y  c  a  represents 

a  qualitative  measure  of  the  forces 

between  the  chains. 

The  authors  would  like  to  thank  Professor  W.  Stockmayer 
for  lending  a  copy  of  Mr.  Hecht’s  thesis  which  gave  a  more 
detailed  account  of  their  work. 
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Apjfflife 

In  Sec,  4,  it  has  been  shown  that  several  new  types 
of  singularities  can  arise  in  idealized  situations  in  which 
valence  forces  resisting  the  bending  of  bond  angles  are  a 
dominant  influence  on  the  spectrum.  For  the  model  considered 
here,  these  arise  from  the  term  x(l  -  cos  cp^)  which  is  approxi¬ 
mately  (xA)^  for  <p^  «  1,  Near  any  critical  frequency  on 
the  face  cp^  =  the  dependence  of  v  on  the  wave  numbers  was 
described  by  an  equation  of  the  form  of  Eq.  (20),  i.e,, 

v2  =  +  aC^  -  <plc)2  +  b(<p2  -  92c)2  +  c(q>3  -  (Al) 

Various  kinds  of  singularities  arose  in  g(v)  depending  upon 
whether  a,  b  and  c  were  zero,  positive  or  negative. 

Table  II  contains  a  list  of  all  the  various  possibilities 
with  c  /  0,  If  c  =  0,  Eq.  (Al)  reduces  to  a  two  dimensional 
quadratic  form  giving  singularities  in  g(v)  of  a  type  already 
considered  by  van  Hove.  Of  the  18  anomalous  singularities 
listed,  many  combine  into  equivalent  pairs  because  the  form  of 
Eq.  (Al)  is  invariant  to  interchange  of  a  and  b  along  with  an 
interchange  of  dummy  indices  1  and  2  on  92  and  92,  thus 
a  <0,  b  >  0  is  equivalent  to  a  >  0,  b  <  0.  Furthermore  if  we 
change  the  sign  of  a,  b  and  c,  this  is  equivalent  to  changing 
v2  -  v2  to  v2  -  v2.  The  18  cases  can  be  grouped  into  pairs 
differing  only  in  that  v  >  v.  for  one  of  the  pairs  is  equivalent 

v 

to  v  <  vc  for  the  other.  A  minimum  is  paired  with  a  maximum, 
for  example. 
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Table  II.  All  possible  types  of  anomalous  singularities  are 
classified  according  to  whether  a,  b  and  c  are  positive  (+), 
negative  (-)  or  zero  (0).  The  fourth  column  indicates  that  the 
singularity  is  a  one,  two  or  three  dimensional  form  (ID,  2D  or 
3D)  of  a  minimum  (min.),  maximum,  (max.)  or  saddle  point  (S.P. ) 
The  fifth  column  indicates  the  approximate  frequency  at  which 
any  of  these  types  appear  in  Fig.  4  (see  Sec,  4). 
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Except  for  these  differences  due  to  interchange  of  9^ 
and  92  or  reflections  of  the  frequency  about  vc,  there  are 
only  six  different  kinds  of  singularities  as  grouped  in  Table 
II,  one  one  dimensional  case,  two  two  dimensional  cases,  and 
three  three  dimensional  cases.  The  present  model  conveniently 
displays,  at  least  approximately,  an  example  of  each  of  the 
six  types,  always  a  case  with  c  >  0  (it  is  hard  to  imagine  how 
c  <0  could  ever  arise  in  a  practical  problem  except  by  accident). 
The  model  also  shows  all  of  the  standard  kinds  of  singularities 
described  by  van  Hove  for  three  dimensional  lattices  and 
approximately  all  the  usual  kinds  of  singularities  for  one  and 
two  dimensional  lattices  as  well. 

Critical  wave  numbers  and  frequencies  of  anomalous 
types  are  to  be  expected  in  quite  general  chain  polymer  crystals 
in  which  the  dominant  forces  resisting  displacements  perpen¬ 
dicular  to  the  chain  result  from  the  bending  of  bond  angles# 

This  does  not,  however,  exhaust  all  the  possible  kinds  of 
anomalous  singularities  of  physical  interest.  A  somewhat 
similar  situation  arises  in  layer  type  crystals  such  as 
graphite^*.  The  forces  resisting  displacements  perpendicular 
to  the  crystal  planes  for  the  lamellar  crystals  seem  to  arise 
also  mainly  from  the  bending  of  bond  angles  but  now  the  fourth 
power  dependence  of  on  wave  numbers  is  associated  with  two 
of  the  three  wave  number  components  instead  of  just  one  as  for 
the  chain  polymer  crystals. 

Corresponding  to  Eq.  (Al),  we  might  obtain  for  the 
lamellar  crystals  critical  points  for  which 
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v2  -  v2  +  a(9x  -  <plc)1+  +  bCcp-L  -  9lc)lf  +  c(93  -  93c)2  (A2) 

Actually  one  is  apt  to  obtain  a  more  general  form 

v2  ~  v2  +  aC©)^  +  “  93c)2  (A3) 

in  which  9  and  9  are  the  polar  representations  of  9-^  -  *LC 
and  92  -  92c> 

91  ~  9lc  -  9  cos  G»  92  “  92c  5  9  sin  ©. 

a(©)  is  some  continuous  function  of  ©  but  not  necessarily  that 
implied  by  Eq.  (A2)  and  9^  can  correspond  to  any  arbitrary 
direction  in  the  wave  number  space  but  in  practical  situations 
would  almost  certainly  be  the  wave  number  component  for  waves 
traveling  perpendicular  to  the  crystal  planes  of  the  lamellar 
structure. 

There  are  not  many  new  types  of  singularities  arising 
from  Eqs.  (A2)  or  (A3)  that  are  of  general  interest.  If  either 
a=0orb=0in  Eq.  (A2),  the  form  of  Eq.  (A2)  reduces  to  a 
special  case  of  Eq.  (Al),  No  physical  situation  comes  to  mind 
that  would  give  negative  values  for  a  or  b  and,  in  fact,  the 
lamellar  crystals  are  expected  to  lead  to  equations  such  as 
Eq#  (A3)  with  a(©)  positive  definite.  There  are,  however, 
still  three  kinds  of  singularities  that  are  certain  to  arise 
in  limiting  cases  of  physically  realistic  situations,  namely 
a(©)  >  0  and  c  >  0  (a  three  dimensional  minimum),  c  <  0  (a 
three  dimensional  saddle  point)  and  c  =  0  (a  two  dimensional 
minimum).  Examples  of  each  of  these  appear  in  the  spectman  of 
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7 

graphite  proposed  by  Komatsu  and  Nagamiya  . 

There  finally  is  the  possibility  that  v2  could  be 

proportional  to  the  fourth  power  of  the  wave  number  components 

in  all  three  directions,  i.e.,  v2  a  homogeneous  fourth  degree 

polynomial  in  the  (9.,  -  9,  ).  It  is  difficult  to  imagine  any 

«»■  j  c 

physical  situation  that  would  give  rise  to  this,  however* 

We  shall  consider  here  only  the  nine  physically  inter¬ 
esting  cases  mentioned  above,  the  six  in  Table  II  and  the  three 
special  cases  of  Eq.  (3A).  Without  loss  of  generality,  we  may 
assume  the  <p2  and  93  are  any  three  independent  wave  vector 
components  scaled  in  such  a  way  that  v2  is  periodic  in  91#  q>2 
and  9^  with  period  2rc.  The  normalized  distribution  density  of 
v2  (for  any  individual  branch  of  the  spectrum)  G(v2)  is  defined 

p  r  p 

G(v2)  =  (2tO“3  dq^dq^dq)^  (A4) 

dv"  JJJ 

integrated  over  the  region 

2  2 
v  (q^j  q>2»  9^)  <  v  ,  0  <  9^  <  2n 

p 

Even  though  the  value  of  G(v  )  does  not,  in  most  cases, 

2 

depend  only  upon  the  local  behavior  of  v  (9^,  92,  9^)  near 

q>j  =  <p^c,  the  singular  properties  of  G(v  )  in  which  we  are 

interested  are  a  consequence  of  the  local  properties  and  we 

may  consider  the  various  cases  of  Table  II  as  if  Eq.  (Al)  were 

valid  for  all  values  of  the  9^'  s.  Either  of  two  things  may 

happen  as  v2 — ^v2.  If  G(v2) — £00  as  v2— >v2,  then  the  value 
2  2  2 

of  G(v  )  as  v — >  is  dominated  by  the  contributions  from  the 
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neighborhood  of  the  critical  wave  vector  where  Eq.  (Al)  is 

2 

accurate.  If  G(vc)  is  finite,  its  value  will,  in  most  cases, 

? 

depend  upon  properties  of  v  (<plf  <p2,  <p^)  far  away  from  the 
critical  wave  vector  where  Eq.  (Al)  is  not  accurate.  If, 
however,  we  evaluate  G(vc)  -  G(vc)  using  the  same  approxima¬ 
tions,  Eq.  (Al),  in  evaluating  each  term,  the  contributions 
far  from  the  critical  point  nearly  cancel  and  we  can  still 
correctly  determine  the  kind  of  singularity  from  the  behavior 
of  G(v  )  -  G(vc)  even  though  we  cannot  find  either  G(v  )  or 
G(v2)  individually. 

p 

From  Eq.  (Al),  we  notice  that  v  is  locally  an  even 
function  of  (<pj  -  9^q)  for  J  =  1,  2,  and  3  so  that  the  con¬ 
tribution  to  G(v2)  from  the  vicinity  of  the  critical  wave 
vector  is  eight  times  that  from  the  region  (q>j  -  <Pjc)  >  0#  If 
we  consider  c  >  0,  which  can  be  done  with  no  loss  of  generality, 
the  range  of  >  0  in  Eq.  (A4)  is 

0  <  <  c^Cv2  -  v2  -  a(<p1  -  <plc)2  -  b(<J>2  -  <P2c ) 2 3^ 

and  so 

0(V2)  -  jj  [v2  -  V2  -  at^  -  9lc)2 

-  b(?2  -  92o)2]'3/4  d<p1d92  (A5) 

The  limits  of  integration  are  automatically  taken  care 
of  by  the  requirement  that  the  integrand  be  real,  i.e.,  the 
quantity  in  brackets  must  be  positive  and  also 

<HV  1  =  1,  2,  3. 


0  <  -  *jo 


(A6) 
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p 

Mj  is  some  number  of  order  1.  If  a  =  0  (v  independent  of  9^) 

p 

we  must  choose  =  it  and  if  b  =  0  (v  independent  of  92)  we 
must  take  Mg  s  Except  for  these  cases,  the  actual  value  of 
Mj  is  unimportant  because  Eqs,  (Al)  and  (A5)  are  only  approxi¬ 
mations  valid  for  19^  -  9-jcJ  <  <  1*  Any  results  that  are 
meaningful  can  not  involve  the  values  of  Mj  for  this  would 
immediately  imply  that  the  results  depended  upon  the  behavior 
of  v2  (91,  <p2,  9^)  in  a  region  where  Eq.  (Al)  is  incorrect. 

The  one  dimensional  result  (a  =  0,  b  =  0)  follows 
directly  from  Eq.  (A5) 


G(  v2)  ~ 


<v2  -  v?)"3/V«lA 


?  2 

for  v  <  v 

c 


_  2  2 
for  v  >  vc 


(A  7) 


with  G(v2)->  oo  as  v2-^  v§. 

For  the  two  dimensional  cases  (b  =  0),  we  make  the 


substitution 


x  =  |a(91  -  91c)2/(v2  -  v|> !1/2 


(A8) 


and  obtain  the  following  limiting  expressions.  For  a  <  0 


(2D  min.), 


G(v2)  =  0 


2  2 
for  v  <  vj 


G(v  )  - 


(v2-v2r*  t 


4n2|a|£c^  (1  -  x2)^1* 


(v2  -  V2)'*  [raA)32 
8x2|2xa|'i'c^ 


P  P 

for  v  >  v* 


(A9a) 


(A9b) 
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Again  G(v2)-->  oo  as  v2~ ^  v2.  For  a  >0  (2D  SP) 


for  v2  <  v2 


2  2 
for  v  >  v 

c 


(A10) 


V  -  vli^cmA)]2  2  , 

- -2 — - -  +  c  for  v*  <  v2 

8n2|7ta|ic^  c 

G(v2)  ~  <  (A] 

|v2  -  v2r*[r(lA>]2  2  2 

- 2 -  +  C  for  v  >  v 

i  l  c 

^  8v?\2m  Pc'* 

in  which  C  is  some  undetermined  constant. 

The  three  dimensional  cases  with  a  and  b  of  the  same 
sign,  (+),  can  be  conveniently  handled  by  a  polar  coordinate 
substitution 

I  a  1^(9^  -  9lc)  =  r  cos  ©,  }bj^(q)2  -  <p2c)  g  r  sin  © 


which  transforms  Eq.  (A5)  into 


G(v2)  -  _ 1 _  _ - 

8it2|ab|^c^  [v2  -  v2  +  r2]^^ 


This  gives  immediately  the  results  for  a,  b  >  0  (3D  min.) 


G(v2)  - 


(v2  -  v2)* 
*f7t2  labile* 


and  for  a,  b  <  0  (3D  SPII) 


G(v2)  -  G( v^)  - 


-(v2  -  v2)* 


for  v2  <  v2 


for  v2  >  v2 
c 


p  p 

for  v  <  vz 


for  v2  >  v2 


(All) 


(A12) 
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Both  Eqs.  (All)  and  (A12)  give  verticle  tangents  on  only  one 
side  of  v2.  ■ 

The  remaining  singularity,  a  >  0,  b  <  0  (3D,  SPIII) 


gives 


G(v2)  -  G(v2)  ~ 
c 


Iv2-  ygl* 
23/2it2|ab|M 

| V2  - 


p  p 

for  v  <  Vq 


for  v2  >  v2 


(A13) 


[,  4it2|  abl^c^ 

which  has  a  (asymmetric)  verticle  tangent  singularity  on  both 


sides  of  v£. 


we  let 


For  the  three  interesting  cases  arising  from  Eq#  (A3), 


(Al4) 


A  =  jL  tC2  d©[a(e)]“1/2 


and  obtain  the  following  expressions.  For  c  =  0  (2D  min.) 


G(v2)  = 


A( v2  -  v2)“1/2 
c 


2  2 

for  v  < 

c 


for  v2  >  v2 


(A15) 


which  has  an  infinite  singularity.  For  c  >  0  (3D  min.) 


G(v2)  = 


|  *“V1/2a 


2  2 

for  v  < 

c 


p  p 

for  v  >  v* 


(A16) 


a  finite  discontinuity.  For  c  <  0  (3D  SP) 


G(v2)  =  C  -  it"1!  cj  ~1/2  A  log |  v2  -  Vq| 


(A17) 


0 


another  infinite  singularity.  C  is  again  an  undetermined 
constant  * 
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gjfittfl.  .gaptloAS, 

Fig.  1.  Three  dimensional  surfaces  of  constant  frequency  In 
the  wave  number  space  (tp-^ ,  <p2>  9^)  for  modes  vibrating  perpen¬ 
dicular  to  the  chains  are  shown  for  several  small  values  of  v. 
For  the  smallest  v,  the  surfaces  are  ellipses  flat  in  the  9^ 
direction.  Due  to  the  valence  forces,  the  surfaces  become 
elongated  in  the  92  direction  for  the  larger  values  of  v, 

Fig.  2.  Three  dimensional  surfaces  of  constant  frequency  for 
larger  values  of  v  are  shown  by  their  intersections  with  the 
pianos  (plv  <p2  and  9^  =  0  or  it.  Solid  lines  are  intersections 
with  the  faces  9^  =  it,  broken  lines  with  the  faces  9^  =  0,  of 
the  cube,  j  =1,  2  or  3*  The  surfaces  are  labeled  with  the 
values  of  v2,  For  small  v2,  the  surfaces  are  elliptical  shaped 
near  the  origin,  see  Fig,  1,  become  elongated  in  the  92  direc- 

p 

tion  as  vc  increases,  become  nearly  cylindrical  for 

2y  <  v  <  a.  For  a  +  4y  <  v2  <  4k,  the  surfaces  are  nearly 

2 

planes  perpendicular  with  the  9^  axis.  The  largest  v  , 

4x  +  a  +  2y,  is  at  the  point  (it,  0,  it). 

Fig.  3.  gx(v),  the  frequency  distribution  for  modes  vibrating 
perpendicular  to  the  chains,  plotted  for  small  values  of  v, 
show  the  transition  from  g(v)«0v2  to  g(v)<Cv^2*  Curve  (1) 
is  proportional  to  v2,  Eq.  (8)}  curve  (2)  is  proportional  to 
v3/2,  Eq,  (9)»  and  curve  3  is  the  accurate  curve  for  g(v), 


Eq.  (10). 
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Fie*  4.  A  plot  of  g^(v)  for  the  complete  range  of  v,  Fig.  3, 
is  a  magnified  view  of  the  small  rectangle  near  the  origin. 

The  solid  curve  shows  an  accurate  curve  of  g^(v).  The  three 
broken  curves,  proportional  to  v^2,  v^2  and  v"1^2  represent 
Eqs.  (9),  (12)  and  (l4),  respectively* 

Fig.  5.  g^(v),  the  spectrum  for  modes  vibrating  parallel  with 

the  chains,  is  plotted  for  the  entire  range  of  v  with  a 
magnified  view  of  the  high  frequency  end  shown  by  the  insert. 
The  broken  curves  are  proportional  to  v2  and  v°  =  const., 

2/n,  corresponding  to  Eqs.  (23)  and  (24),  respectively. 

Fig.  6.  The  complete  frequency  distribution  for  all  modes, 
g(v),  is  shown  by  the  solid  curve.  The  step  function  curve 
is  the  approximate  curve  computed  by  Stockmayer  and  Hecht 
for  the  same  system. 
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